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1. Introduction 

Intrinsic localized modes (ILMs or discrete breathers) have been a center of intense 
theoretical, numerical, as well as experimental investigations over the past two decades; 
see e.g. the reviews [U El [3l Hj. Since their theoretical inception in the context of 
anharmonic nonlinear lattices [5j E] and subsequent rigorous proof of existence, under 
appropriate nonresonance conditions [7], numerous experimental realizations of such 
structures have arisen in settings ranging from optical waveguides and photorefractive 
crystals, to micromechanical cantilever arrays and superconducting Josephson junctions, 
as well as Bose-Einstein condensates and electrical lattices, among many others [3]. 

While the most fundamental modes among these discrete breathers, namely the 1- 
and 2-site solutions have been analyzed in some detail, much less is known about the 
case of multi-site breathers or multibreathers. The latter were initially discussed in [7j. 
Since that pioneering work, a lot of effort has been invested in proving the existence 
of multibreathers in Klein-Gordon chains (e.g. [El El HD])- Regarding the stability of 
these motions, in [11] some stability results are obtained using the formulation of [10J, 
but these results are applicable only to few-site excitations. On the other hand, in 
[T2] (see also the more recent discussion of [13]) some theorems about the stability of 
multibreathers are proven, using Aubry's band theory [I], which can be applied to an 
arbitrary number of site excitations. 

In the present work, a n-site breather stability theorem is proven. It generalizes 
the two previously mentioned works, by proving a detailed counting result about the 
number of real and imaginary characteristic exponents of the corresponding breather for 
arbitrary configurations. It should be noted that this result proves a relevant statement 
made in [12] as a claim based in numerical findings. The relevant eigenvalues are 
estimated to 0{\fe). Our method is based on the notion of the effective Hamiltonian 
originally introduced in [H] and generalized in [151 [IB]- This idea has already been 
used in order to prove existence and stability of multi-site breathers in hexagonal and 
honeycomb lattices [171 03 H2!- Similar results have been acquired for the case of the 
discrete nonlinear Schrodinger (DNLS) lattice [201 [22] and were recently used in the 
study of discrete solitons in hexagonal and honeycomb lattices in [23] . 

Our principal result shows that for soft nonlinearities and general, multi-site 
excitations, the relevant structures may only be stable (for positive values of the 
coupling) when the adjacent sites are out of phase by n with respect to each other. 
For negative values of the coupling, stability is possible for in-phase excitations. This 
situation is reversed in the case of hard nonlinearities (i.e., in-phase multibreathers are 
stable for positive weak coupling, while out-of-phase ones for negative, weak coupling). 

Our presentation is organized as follows: in section [2] we define the system under 
consideration, in section [3] we set up the general conditions for existence of multibreather 
solutions, in section 0] we acquire the theorem about the stability of the previously 
mentioned solutions and finally in section [5] we perform some numerical calculations in 
order to verify our theoretical predictions. 
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2. Definition of the system - Terminology 

We define our oscillators by an autonomous Hamiltonian of one degree of freedom 
H u = l -p 2 + V(x) 

where V(x) is the potential function. In this case, the system is integrable since H u 
is always an integral of motion. We assume that V(x) possesses a minimum at x = 
(without loss of generality) with V"(0) = uj^ with u p G R. 

Because of the time reversal symmetry x(—t) = x(t),p(—t) = —p(t) the solution of 
the oscillator can be written as 

oo 

cos(nw) (1) 

n=0 

where J, w are the action-angle variables. Note that in the action-angle variables the 
motion of the oscillator is described by 

w(t) = UJt + w 
J it) = const. 

where to is the frequency and Wq is the initial phase of the periodic motion. 

We construct our chain by considering a countable set of oscillators with a nearest- 
neighbor coupling through a coupling constant e. The Hamiltonian then becomes 

OO / * \ oo 

H = H + eH l = Up" + + \ E ^ - ( 2 ) 

i=— oo ^ ' i=— oo 

where Xi is the displacement from the equilibrium, and pi the momentum of the i-th 
oscillator. Note that Hq is trivially integrable, being separable. 

3. Existence of multibreathers 

Consider the "anticontinuous" limit e = where n + 1 adjacent "central" oscillators 
move in periodic orbits with frequency u but arbitrary phases, while the the remaining 
"non-central" oscillators lie at rest (xi,pi) = (0, 0). This state defines a trivially localized 
and time-periodic motion with period T = 2-k/uj. We seek conditions under which this 
motion can be continued for e ^ to provide a multibreather of the same frequency u. 
In the next section we will determine the linear stability of the resulting solutions. 

We apply the action-angle canonical transformation to the central oscillators. The 
system is described now by the set of variables (xi,pi, Wk, Jk) with k G § and i G Z \ S 
where S is the set of "central" oscillators. So the periodic orbit which corresponds 
to the multibreather is described at time t by z{t) = (xi(t),pi(t),Wk{t), h{t)) with 

z(t + T) = (x i (t), Pi (t),W k (t) + 27T,I k (t)). 

In [H] (extended in [15] and [16J) it is proven that under the non- resonance 
condition nu ^ u^in G Z there is an effective Hamiltonian H eG whose critical points 
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correspond to periodic orbits (in fact, breathers) of the full system for e small enough. 
The effective Hamiltonian is defined by 

# eff (I;,A<&) = ^j>Hoz{t) dt, 

where z is a periodic path in the phase space obtained by a continuation procedure for 
given relative phases (pi, relative momenta I and symplectic "area" A. In the lowest 
order of approximation, the unperturbed orbit z can be used instead of z. In our case, 
this coincides with the averaged Hamiltonian over an angle, for example wq = ut + wo , 
due to the linear relationship of wo with t. Since, by construction, the resulting effective 
Hamiltonian does not depend on the selected angle w , and due to the nature of the 
system, a canonical transformation to the "central" oscillators is induced 

$ = u>o A = J Q + . . . + J n 

4>l = Wi - w h = Ji + . . . + J n 

(p 2 = W 2 -W 1 I 2 = J 2 + ■ ■ ■ + Jn (3) 

4>n = W n - U> n _i I n = J n 

and, in the lowest order of approximation, the effective Hamiltonian becomes 

H eS = H (I i ) + e(H 1 )(<f> i ,I i ) i = l...n (4) 

with 

(H 1 ) = ^H 1 dt 

where the integration is performed along the unperturbed periodic orbit. Note that 
(Hi) coincides with (Hi) WQ , the average value of Hi over the angle w and since H eS is 
independent of A is a constant of motion. 

As we have already mentioned, the critical points of this effective Hamiltonian 
correspond to breathers. But for non-degenerate critical points, to leading order in e 
this condition reduces to the conditions (which were aquired also in [ID] ) 



d(Hi) _ d>(Hi) 



d 2 Hn 



dJidJj 



7^0^ |^0, u p ^ku. (5) 



Note that, since we consider central oscillators with the same frequency, condition (J5b) 
can be reduced to |j ^ 0. By taking into account ([I]) we get (see Appendix A| for 
details) 

1 oo n 

^> = -2EE^ cos w*) ( 6 ) 

m=l i=l 

hence, the condition (jSK) becomes 

oo 

2j mA 2 m sin(m0j) = (7) 



m=l 
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which has at least the solutions 

(pi = 0,7T. 

Intuitively, we expect that these are the only solutions, in this kind of systems, however, 
a general proof of this conjecture is not presently available [21]. 

It is interesting to note in passing here the similarity of the above conditions to the 
Lyapunov-Schmidt persistence conditions obtained in the context of the DNLS model 
in [201 [2B [22]- However, in the latter case, the presence of a single frequency enforces 
the (pi = 0,ii condition. 

Remark 1. Although we don't have a full proof of the above assumption yet, physical 
considerations suggest its potential validity. For instance, consider the simplified setting 
wherein the displacement x(w) (equivalently x(t)) from equilibrium is described by 
the truncated series x(w) = A + A\ cos(w) + A 2 cos(2u>). Then, the acceleration 
a(w) = x{w) reads a{w) = — u 2 [A\ cos(w) + 4A 2 cos(2w)]. Since we know that a 
in the two edges of the motion should be a(0) < and aiit) > 0, this means that 
A 1 + AA 2 > and A 1 - AA 2 > and finally A\ > 16A 2 .. Since for this case (C]) reads 
A\ sin(0) + A\ sin(20) = sin(0)[yl^ + 4v42cos(0)] = 0, we conclude that in this special 
case, the above physical considerations preclude solutions other than <p = 0, it. 

4. Stability of the mult ibr eat her solutions 

The linear stability of the fixed point of H determines also the linear stability of the 
breather. This is proven in [H] for the first order approximation to H cS , under the 
assumption of distinct eigenvalues of the first order matrix, and in [15] for the general 
case. This fact has already been used in order to study the stability of 3-site breathers 
in [18]. Again, there is a direct analog of this in the DNLS case, whereby the Jacobian of 
the Lyapunov-Schmidt conditions in [2D1 EH [22] is, to leading order, directly analogous 
to the squared eigenvalues of the full linearization problem. 

To make things more precise, the linear stability of a multibreather, is determined 
by its Floquet multipliers (see e.g. [1]), which are the eigenvalues of the monodromy 
matrix of the corresponding periodic orbit. If all the multipliers lie on the unit circle the 
breather is linearly stable, otherwise the breather is unstable. Due to the Hamiltonian 
character of the system if A is a multiplier so are A*, A -1 , A* - . In particular, for 
multibreathers, when e = all the multipliers lie in two conjugate bundles at e ±luJpTb 
except for n + 1 pairs, which lie at unity and correspond to the central oscillators. When 
this solution is continued for e ^ the multipliers which belong to the two bundles, 
being of the same Krein kind (e.g. [22]), move along the unit circle to form the phonon 
band. On the other hand, one pair of the multipliers of the central oscillators will 
remain at 1 because of the corresponding invariance of the system while the rest can 
move either along the unit circle or outside the unit circle determining in this way the 
linear stability of the multibreather. 
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We define the characteristic exponents o~i of the multibreather, or equivalently of 
the corresponding periodic orbit, as 



A; 



The non-zero characteristic exponents of the central oscillators correspond to the 
eigenvalues of the (2n x 2n) stability matrix [HI [15] E = flD 2 H cS , where fl = 

-I 

1 O 

we demand that all the eigenvalues of E be purely imaginary. The stability matrix E, 
to leading order of approximation and by taking into consideration becomes 



and I the nxn identity matrix. According to the above, for linear stability 



/ 



E 



A 


B 


C 


D 



eAi 


eB 1 


Co + eCx 





d&idh 



d 2 H { 
\ dljj 



dltdlj 



Using (jSJ), the elements of E take the form 



— — = ^m 9 (J)sm 



m<t>i 



m=l 







with g{J) = — (A m (Ji-i)A m (Ji)) 
while, (|Appendix B[), 



--Ji=J 



and 



hJ 



d 2 m \ 



'jdli J 



n. 



(8) 



d 2 H 
dhdh 



i o duJ 




a J 


3 = 


< duo 




~dJ 


3 = 





else 



and 



<9 2 (#i) 



fib 




3 = i 
3 ^ i 



where 



/(0) = ^E n2A « cos ^)- 



(9) 



Since we consider solutions with 
dHl) becomes 



0, 7r, where A = D = O, the stability matrix of 



E 



O B 

C O 

d 2 (H x ) 
dhdh 



O ^Bi 
Co + ^Cx O 



(10) 



terms contribute only to higher (than the leading) order, 



Due to Lemma [TJ the 
and hence will not be considered further in what follows. 

Lemma 1. The leading order of approximation of the eigenvalues o/E is O(yfe). The 
term Ci = 9 d j^ in fifty only affects the eigenvalues at 0(e 3 ^ 2 ). 
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As it is also shown in Appendix C up to the leading order of approximation, we have 
a ±i = ±^xTi + 0(£ 3/2 ) i = ] 



n. 



(11) 



where a±i (the characteristic exponents) are the eigenvalues of E and xu are the 
eigenvalues of BxCV So, the sign of Xu define the stability of the multibreather. Let, 
fa — Then, BxCo becomes 

/ Vi -h o 

-/a 2/ 2 -/ a 



Bi ■ C 



duj , 
dJ' 



du 
"dl 



\ 











V 



(12) 



fn—l ^fn—1 fn—1 

Lemma 2. Let Zj 6e tae eigenvalues ofZ. Then, the number of positive Zi's equals the 
number of positive fa's, while the number of negative Zi's equals the number of negative 
fa's- 



The proof can be found in |Appendix D 



Lemma 3. Assuming the absence of solutions of |3]a) other than < 
and f(n) < 0. 



□ 

0,7T ; thenf{0) > 



Proof. The fact that /(0) > is obvious from ([9]) since A{ are the Fourier coefficients 
of a smooth real function. On the other hand, 



F(4>) 



cos nc 



n=l 



is a continuous function. Since the values = 0, rr correspond to the extrema of 
F((j)), because of continuity, one of them corresponds to a local minimum while the 



other corresponds to a local maximum. So, since /(0) = d A > corresponds 

4>=o 

to the minimum, therefore /(</>) = vr must correspond to the maximum of -F(</>) and 

f M = q» <o. 



Lemma 4. I/efy < and . 



7i Vi 



n, or i/etr > and . 



iaen a// iae eigenvalues of E are purely imaginary up to 0(\fe) terms. 



□ 



Vi = 1 . . . n, 



Proof. Due to ffi2|) . we have x 



— So, by using ffTTl) we get 



±\/- £^^ + C(e 3/2 ) 



(13) 



The sign of Zi is defined by the value of <j>i according to lemmas [2] and El which, in turn, 
completes the proof of the lemma. □ 
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If the eigenvalues Aj are imaginary and distinct up to 0{y/e) terms the higher order 
terms cannot push them outside the imaginary axis for a variance of e, say Ae, small 
enough, because of continuity. On the other hand, if the eigenvalues Aj have multiplicity 
> 1 up to 0(-\/e) terms the higher order terms can, in principle, push them outside the 
imaginary axis for Ae arbitrary small, which would cause complex instability, through 
a Hamiltonian Hopf bifurcation. This, however, cannot happen in our system since a 
specific symplectic signature property holds. 

Lemma 5. If the eigenvalues of E are imaginary up to 0(y/s) terms they remain 
imaginary up to all orders of approximation. 

Proof. If the eigenvalues of E are imaginary up to some order of approximation, then, 
according to [26], if the corresponding quadratic form of D 2 H eS is definite, then the 
eigenvalues remain imaginary for all orders of approximation. The matrix D 2 H cS is 



/ d 2 H Q \ 



D 2 RC s 



OJidJj 
O 



and the corresponding quadratic form is S 2 H eS = (I, 0) ■ D 2 H cS ■ (I, 0) T , with / = 
(Ji, . . . , J n ) and = . . . , 4> n ). Finally we get 

^ = fj ^ + ih ~ /l)2 + ' ' ' + ( J « " J -l) 2 + ^ + 8 IffaM + --- + / • 

This quadratic form remains definite for all the configurations which are described in 
Lemma HI So, even in the case of higher multiplicity, the imaginary eigenvalues of E 
remain on the imaginary axis. □ 

The sequence of the above lemmas leads to our main stability theorem, as follows: 

Theorem 1. Under the assumption that |?|) has no other solutions than 0« = 0, n, then, 
if Ejjj < the only configuration which leads to linearly stable multibreathers, for |e| 
small enough, is the one with <fii = ir Vi = l...n (out-of-phase multibreather) , while 
if ETjj > the only linearly stable configuration, for \e\ small enough, is the one with 
4>i = Vi = l...n (in-phase multibreather). Moreover, for e^j < (respectively, 
e|j > 0), for unstable configurations, their number of unstable eigenvalues will be 
precisely equal to the number of nearest neighbors which are in (respectively, out of) 
phase between them. 

Proof. Since for a linear stable multibreather we need imaginary eigenvalues of E, the 
only possible configurations for stability are the ones described by the theorem, as it can 
be shown from lemmas [3] and HI The multibreather will remain stable for small enough 
values of |e|, until the eigenvalues which correspond to the central oscillators will collide 
with the linear spectrum, causing a Hamiltonian Hopf bifurcation, leading to complex 
instability. □ 
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We note in passing that the above theorem bears a direct analogy to Theorem 3.6 
of [20] for the DNLS case. 

Remark 2. Note that if the the on-site potential is even V(x) = V(—x) then the cosine 
Fourier series of x(t) becomes 

oo 

x(t) = A + Y^ Mn-i cos[(2n - l)w] 

n=l 

and /(0) becomes 

f(<f>) = 2 E( 2n - l T A ln-i ™s[{2n - 1)0] 

n=l 

which means f(ir) = — \ J2^=i(2n ~ l) 2 ^L-i < 0- So, Theorem [T] can be reformulated 
without the need of exclusion of possible other solutions of (JTj). 

5. Numerical Results 

As a prototypical numerical demonstration, consider a chain consisting of oscillators with 
on-site quartic potentially) = ^—0.27^—0.03^. This potential is softening (|j < 0) 
as it can be seen in figure [TJ We will consider the orbit with period T = — = 7.434 
which corresponds to amplitude of oscillation x max = 1.949275 =^ J = 1.20306 =>■ |j = 
-0.224556. For the same orbit we get /(0) = 1.423404 and /(tt) = -1.279544. 




(a) (b) 

Figure 1. Panel (b) shows the dependence of u> with respect to J for the specific V{x) 
of panel (a). 



5.1. 2- site multibreathers 

We consider first the case of two cental oscillators (two oscillators moving at the anti- 
continuous limit). In this case there is only one <ft = w 2 ~ Wi — ^20 — ^io an d 
consequently only a pair of characteristic exponents. The leading order approximation 
of o~i is, according to (jTBl . 



o ±l = ±\l2e d ^f{<P). (14) 
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The resulting 2-breathers are: 

• The in-phase 2-breather with = 0. In our example where e > and |j < 
we get er±i £ R, which leads to an unstable 2-breather. In figure Eh, the profile 
of the in-phase 2-breather is shown while in figure [2b the real part of the positive 
characteristic exponent of the central oscillator o\ as calculated by the numerical 
simulation is shown (solid line) together with the theoretical 0(y/e) prediction of 
o~i in ([14"]) (dashed line). We can see that for small values of e the agreement is 
excellent, while for larger values of e, where the higher order terms of o~i become 
significant, the two lines start to diverge. 

• The out of phase 2-breather with = it (see figure [3]). In our example it is o±\ G I, 
which leads to a linearly stable 2-breather. In figure EH, the profile of the out of 
phase 2-breather is shown, while in figure [2b the imaginary part of o\ is shown. 
The solid line represents the numerically calculated value while the dashed line is 
the theoretically approximated value. For small values of e the two lines again 
nearly coincide, while for larger values of e, the two lines diverge. For e ~ 0.0254 
the solid line possesses a cusp which results from the collision of o\ with the linear 
spectrum. At this point two characteristic exponents (and their conjugates) acquire 
a nonzero real part, through a Hamiltonian Hopf bifurcation, and the corresponding 
multibreather becomes unstable. This is the typical mechanism through which 
multibreather solutions identified herein as stable for small e eventually become 
unstable as the coupling is increased. 




- 10 



10 



X) 



0.C05 



0.010 



0.015 



0.020 



0.025 



0.030 



sites 




(b) 



Figure 2. (a) The profile of an in-phase 2-site breather for e = 0.02. In (b) the real 
part of a\, for increasing values of e, is shown. The solid line represents the numerically 
calculated value, while the dashed one is the one resulting from (|14p . 



5.2. 3 -site multibreathers 



The next step is to consider three central oscillators. In this case there exist two 
independent 0j's, <f>\ and 02- So, there are three relevant configurations to examine, 
which correspond to the three possible combinations of 0j. 
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Figure 3. (a) The profile of an out of phase 2-site breather for e = 0.02. In (b) the 
imaginary part of a±, for increasing values of e, is shown. The solid line represents the 
numerically calculated value, while the dashed one is the one resulting from (fl4|) . 



0i — 02 = (In-phase multibreather) . Following (fTBl . the leading order 
approximation of the four characteristic exponents of the 3-breather is 



^ = ± A /- e ^/(0) 



cr ±2 = ± 



-3e>, 



(15) 



In our example this is an unstable configuration since <J±\,±2 £ In figure UK the 
profile of a 3-site in-phase breather for e = 0.02 is shown, while in figure Hb the 
real part of the corresponding characteristic exponents eri^ is shown. Again, the 
solid line denotes the numerically calculated values and the dashed ones represent 
the theoretical predictions. The agreement is very good, especially for small values 
of e, illustrating the accuracy of our theoretical predictions. 






sites 



(a) 




0.015 0.020 

E 



(b) 



Figure 4. (a) The profile of an in-phase 3-site breather for e = 0.02. In (b) the 
real part of ai,2, for increasing values of e, is shown. The solid lines represent the 
numerically calculated values, while the dashed lines correspond to the theoretical 
result of (fl5l). 



0i = 02 = vr (Out-of-phase multibreather). In this case, the leading order 
approximation of the corresponding characteristic exponents is 

*±i = ±yj-e^f{-K) , a ±2 = ±yl-3e^f(n). (16) 
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In our example this is a stable configuration since cr±i,cr±2 G I. In figure the 
profile of a 3-site out of phase breather is shown, while in figure Hb, the imaginary 
part of the corresponding characteristic exponents is shown. Once again, the 
agreement between the two lines can be noted, at least for small e. For e ~ 0.019, 
the line which corresponds to o~i appears to change slope, a feature which is due to 
its collision with the multipliers stemming from the phonon band. For larger values 
of e the multibreather is unstable. 




(a) (b) 

Figure 5. (a) The profile of an out-of-phase 3-site breather for e — 0.02. In (b) the 
imaginary part of CT1.2, for increasing values of e, is shown. The solid lines represent 
the numerically calculated values, while the dashed ones are the ones resulting from 

(USD. 



• 0i = 0, <p2 = 7T- In this case, the corresponding leading order approximation of the 
characteristic exponents is 

*±i = ±\/-ef7 Ui + 7* - Vf? + fi - JJ^) 

1 — i (17) 

o-±2 = ±a/-£§ [fx + f2 + VA 2 + fl ~ 7J2) 

with fi = /(0) and = f(jt). In our example, a±± G I and a±2 E R [it is 
straightforward to show that this will always be the case if fif 2 < 0], so the 
corresponding configuration, shown in figure El is unstable. In figures [6(b) and 
HKc), the imaginary and real parts of o\ and o~2 are shown. Note that for e ~ 0.03, 
<7i enters the phonon band as can be seen in figure [6](c), where its real part becomes 
nonzero and the corresponding multibreather becomes unstable. 

5.3. 5-site multibreathers 

Our methodology can be numerically applied (through the simple numerical calculation 
of the eigenvalues of a n x n matrix), even when we cannot analytically calculate the 
eigenvalues of E (Z). In this case we can still numerically calculate the eigenvalues of 
Z and get the 0(y/e) prediction from ffTBl . In order to demonstrate this, we consider 
five central sites, so there exist four independent (relative angles) <fi. A representative 
configuration is = <p 2 = and 3 = 4 = n which results in an unstable multibreather 
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(a) (b) (c) 

Figure 6. Panel (a) shows the profile of the multibrcathcr corresponding to the 
4>i = 0, <p2 = 7r configuration for e = 0.02. In (b) and (c) the imaginary and the 
real parts of eri,2, for increasing values of e, are shown correspondingly. The solid 
lines represent the numerically calculated values, while the dashed ones are the ones 
resulting from (fT7|) . 



with Oi 2 G M and 03.4 G I, as expected from our main theorem above (see e.g. figure 

ED- 




Figure 7. In (a) the profile of the 5-breather corresponding to the 4>i,2 = 0, 03, 4 = it 
configuration is shown for e = 0.02. In (b) and (c) the imaginary and the real parts 
of <7j, for increasing values of e, are shown, respectively. The solid lines represent 
the numerically calculated values, while the dashed ones stem from the numerical 
calculation of the eigenvalues of the matrix Z. 



6. Conclusions 

In the present paper, we proved a linear stability criterion for ra-site multibreathers. 
This result generalizes the ones acquired in [11], which can determine the stability only 
for configurations up to three "central" oscillators, while it proves a counting result 
concerning the number of real and imaginary characteristic exponents of the breather, 
which is similar to the claim stated in [12] . In addition, our approach provides an 
O(tJe) estimate of the characteristic exponents of the multibreather solution. Finally, 
the numerical simulations showed that our estimate is accurate for small values of e, 
while it diverges for larger values of the coupling constant, which is naturally expected, 
since for this range of values the higher order terms of the expansion of the characteristic 
exponents become significant. 
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It would be especially interesting to extend these considerations to higher 
dimensional settings, to obtain a systematic characterization of vortex solutions and 
their stability, in square, as well as non-square geometries. Such efforts are currently in 
progress [T7] and will be reported in future publications. 



Appendix A. Calculation of (Hi) 

The average value of Hi is defined as 
(Hi) = ± <f Hidt, 

where the integration is performed along the unperturbed periodic orbit. Since in the 
anti-continuous limit 6 = the only moving oscillators are the "central" ones, Hi 
becomes 

n n 

Hi = x\ - yxjXj-i. 

i=0 i=l 

Therefore, only the mixed terms of Hi interest us since, as we can easily conclude 
following the procedure below, the integration of the square terms over a period provides 
constant terms, i.e. terms independent of fa. We define U = f Q XiXi-idt, so 

1 n 

(Hi) = —— ^J^i + (independent of fa terms). (A.l) 
i=i 

Since at the anti-continuous limit the motion of the oscillators can be described by ([1]) 
we get, by dropping the constant term of the fourier series, for Ii 

h = Xix dt = / \^\^A m (Ji)cos(mwi)A s (J )cos(sw )dt = 

J ° J ° m=ls=l 

oo oo 

= A m ( Ji)A s ( J ) / cos[m(ut + Wi )] cos[s(ut + w 00 )]dt = 

m=l s=l ^° 

~ A m A s f f T 
= 7^ 7^ ™ < / cos[(m + s)ut + (mwio + sw 0Q )]dt 

m=l s=l 2 

+ J cos[(m — s)ut + (mwi — sw 00 )]dt^ . 

Without loss of generality, we can impose that m,s > 0. Then, the only terms that 
survive are the ones with m = s, so we get 

h = -£r I cos[mOi -woo)]d£ = ^ / cos[m(^i-w )]dt = ^ — ^cosmfa. 

m=l ^° m=l ^° m=l 

So, by I 1A.1I) we get 

. n oo 



i=l m=l 
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dH 

Appendix B. calculation of the — — terms 

By taking the inverse transformation of ([3]) we get for the Jj's 

Jo = A-h 
J i = h-h 
J2 = h — h 



(B.l) 



Jn—1 Ifi—1 



Since the integrable part of the Hamiltonian is written by definition as H 
H (Jq . . . J n ) and consequently H — H (Ii . . . I n ) we get by using (IB. II) 



dH _ dH dJ^ dHodJ, dH dH 

~di~ ~ dj~ 1 ~di~ + ~dj~w t ~ ^dj-; + ~dj~- + ^w, 

where the last equality holds because H is separable being the sum of one degree of 
freedom Hamiltonians. So, every frequency depends only on the corresponding action. 
So, we have 

d 2 H d dH d dH _ d 2 H dJ^ d 2 H dJj _ d 2 H d 2 H 
81? JfiidJi-x + dTi~dJ~~ dJ}_ x dh + dJf~dTi~ dJf_ x + DJf 

and for Ui = u we get 

d 2 H 2 du 



dl 2 ~dJ' 
Using the same arguments we get 

d 2 H d dH d dH d 2 H dJi d 2 H 



dl i+1 dli dl i+ i dJi-i dl i+ i dJi dJf dl i+1 dJ\ 
which can be written as 

d 2 H _8lj_ 
dl i+1 dli dJ' 

Appendix C. Expansion of the eigenvalues of E 

Let 

E = 






B 


c 


O 



be the stability matrix, with eigenvalues ov The eigenvalue problem for this matrix can 
be rewritten as 

iBC-a 2 !! =0 
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or 

|BC- X I| = 0. (C.l) 
Using the expansions Xi — Xoi + £Xu + £<2 X2i, B = eBi and C = Co + eCi, we get 

I^BxCo + e 2 BiCi - (xo + e X i + e 2 X2 )I\ = 0. 

Since this relation must hold in the limit e — > 0, we get 

|XoI| = Xoi = i = l... n. 

Hence, condition (1C.1I) becomes 

I^BxCo + e'Bid - (exi + ^X 2 )I| = 

or 

| £ [B 1 C o -XiI + £(B 1 C 1 -x 2 I)]| = 

or 

e n iBiCo - Xil + ^(BiCx - X2 I)| = 0. 

For £ 0, this becomes 

IBjCo - Xil + e(BiCi - ux2l)\ = 0. (C.2) 
However, once again, this condition must hold for the e — > limit which reads 

|BiC -XiI| = 0. 

So, xii depend only on Bi and Co, while Ci only affects the higher order terms. Since 

o- 2 = £Xi + £2 X2 



V Xi 

a = ±V^Xl(l + £*2 + . . .) 

2xi 

hence, up to terms 0(e), the eigenvalues of E are determined by B x and C , while the 
influence of Ci moves to terms of 0(e 3 ^ 2 ). 

Appendix D. Positive and negative eigenvalues of Z 

For reasons of completeness, we also present a proof of the fact that the number of 
positive eigenvalues of Z (zj) equals the number of positive fiS and the number of 
negative Z{ equals the number of negative /jS. This can be done by induction, directly 
following the steps of Appendix C of [27J. 
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First we define Z„ as 



/ Vi -fx 
-h 2/ 2 







V 



-h o 

'fn—l ^fn—l fn—l 
"/n 2/ n / 



(D.l) 



The determinant of Z n is given [27] by 



detZ„ = det 



/ 2/i 

— /2 |/2 



V 



n-1 



n ^ 
l -lJn- 



(n 



+ l)f[/< 



f »+i f / 

Jn n Jn / 



The claim holds for Z, with i = 1, 2, 3. Let's assume that it holds for Z n _x- We will 
examine if it holds for Z n . Note that, we consider only the case fi 7^ 0, since in order 
to have fi = 0, special symmetry conditions should hold. 

Let's consider /j 7^ for i — 1 . . .n — 1. We define f — (fi - ■ ■ f n -i) and / = (/, e) 
with eel, then the eigenvalues Z\(e) . . . z n (e) of Z n are C 1 in e; see e.g. [25] . 

Consider e = first. Since fi 7^ for i = l,...,n — 1 we have that 
Zi(0), . . . , z n _i(0) 7^ 0. In addition 



n-l 



n-l 



'n-l 



'n-l 



Y[ Zl (0) = nl[ ^ ^ sign I Y[ Zi (0) ) -sign JJ/J ^0 



i=l 



i=l 



,i=l 



,i=l 



and £ n (0) = since the last row of Z n vanishes. For e 7^ we have 



n 



n-l 



Zi(e) 



[n 



so, for small e it is 



But, 



sign(> n (e)) = sign(e). 



n-l 



detZ n = (n + l)eJJ/ i ^0 



is valid for every e 7^ 0, so no eigenvalue can change sign as long as e is nonzero. Con- 
sequently the claim holds for Z n which coincides with Z, so the lemma is proven. 
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